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$_i Abstract 
-t— > 

CZ) We consider the XXZ quantum spin chain in its massless, disordered regime at zero field. We 

derive an exact expression for the two-spinon form- factor of S z — \a z by taking a limit of the 
massive XYZ form-factors found by Lashkevich and by Lukyanov and Terras. This result is used to 
find the two-spinon contribution to the spectral decomposition of the longitudinal structure factor 
S zz {k,w). We find that this contribution provides an accurate approximation to the full structure 
factor over a wide range of the anisotropy parameter. The asymptotic behaviour of S zz (k,w) is 
computed as the upper and lower w thresholds of the two-spinon (w, k) band are approached, and an 
analysis of the region of validity of this threshold behaviour is performed. Our results reproduce and 
t-H refine existing threshold behaviour predictions and extend these results to an accurate description 

throughout the two-spinon continuum. 

Dedicated to Professor Michio Jimbo on his sixtieth birthday 

1 Introduction 

1— ( Interacting quantum systems have unique properties when space is one-dimensional (TJ . On the one 
hand, the simple fact that particles cannot avoid each other means that the nature of quantum 
dynamics is complicated in one dimension. The inevitability of particle interactions means that all 
excitations are collective ones. In particular, the quasi-free excitations of Fermi liquids in higher 
dimensions are replaced by the non-perturbative excitations described at low energies by the theory 
of Tomonaga-Luttinger liquids in one dimension. On the other hand, the tools for dealing with non- 
perturbative systems are far more sophisticated for one dimension than for higher dimensions. For 
the class of systems that are quantum integrable, the mathematical toolbox is particularly full. 

The most studied interacting, one-dimensional quantum integrable system is the Heisenberg, or 
XXZ, quantum spin chain j2j 13] 

I N 

H XXZ = - i + oMw - Aofofo) (1.1) 
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with A = cos(7r/(£ + 1)). This model has two nice properties: many exact non-perturbative results 
exist for both finite and infinite N (see for example |H E] and the many references they contain) ; and 
it is experimentally realised. When J > 0, the XXZ chain has a massive antiferromagnetic phase 
for A > 1 and is realised for example by CSC0CI3 [6J. When |A| < 1 the model has a massless 
disordered phase and has been realised experimentally by frustrated spin ladder systems El [9] and, 
very recently, has become in principle accessible using optical lattices |10j-|14j. 



In paper [T5] , we have considered the N — > 00 limit of the Hamiltonian (1.1) in the massless phase 
with < A < 1 (the change in gauge from Hxxz to the Hamiltonian H considered in [15] is given 
in Section 4 of the current paper). In this phase, the system is a Tomonaga-Luttinger liquid |X6|, I17j 
whose fundamental particles are 'spinons': spin-1/2 excitations that can be viewed as domain walls 
dressed with quantum fluctuations [18]. When A = these excitations are non-interacting and are 
described by free fermions. Away from A = the spinons are shaped by the interactions in the bulk, 
and these interactions can be probed by determining how spinons contribute to correlation functions. 
The correlation function we have considered in detail in [TS] is the longitudinal structure factor (LSF) 

/oo 
dte ia; *(vac|5|(i)^(0)|vac) (1.2) 

where S z = \a z . The LSF can be measured directly in neutron scattering experiments (see [19] and 



2 l 

references therein). S zz (k,ui) can be computed by resolving the identity in terms of a complete set of 



spinon states I = ^ an d inserting into (1.2) to give the spectral decomposition 



S zz (k,w) = ^(2tt) 2 <5(£; - K(a))S(w - W(a) \( vac| <Sq I ct) I , 

a 

where K(a) and W(a) are the momentum and energy of the state \a). In our paper |15j . we have 
presented the result for the exact two-spinon contribution to this sum and shown that this contribution 
is a highly accurate truncation, saturating two independent sum rules to around 99% at A = 0.5. 

The main purpose of the current paper is to explain the derivation of the results presented in brief 
in |15j . In particular, we show how the relevant two-spinon form factors (vac|Sg|a!) are obtained for 
the massless phase of the XXZ model. The general method that we follow to obtain massless XXZ 
form- factors is usually called the vertex operator approach (VOA) . The VOA for the antiferromagnetic 
XXZ model is described in detail in [20] , where the representation theory of the quantum affine algebra 
Uqisl-z) plays an essential role. This theoretical framework has been exploited to offer results on 
dynamical correlation functions of the Heisenberg chain both at the isotropic antiferromagnetic point, 
where two |21| [22] and four-spinon [21"] contributions have been obtained, and for the gapped 
antiferromagnet, where two-spinon contributions to the transverse correlator were given [25, 26j. In 
order to extend this approach to deal with the massless regime we follow the strategy proposed in 
[27J: we use the VOA for the XYZ spin chain in the principal regime, map it to the XYZ disordered 
regime, and then take a massless limit to the XXZ model. 

The XYZ model Hamiltonian is given by 

The VOA for the XYZ model was developed in [281 [29]. The role of U q (si 2 ) in the XXZ model is taken 
in the more general XYZ case by the elliptic quantum group of vertex type A q<p (sl2) [301 EH E2]. 
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The VOA is directly valid in the principal regime of the XYZ model for which \J y \ < J x < —J z . 
However, as demonstrated in [33], is is possible to map the principal regime to any other region in 
the phase diagram of the XYZ model. In particular, we can map to the disordered region for which 
I Jz I < Jy < Jx • A transformation which achieves this is H^yz d = ^ ^xyz ^ > wnere 



G = ■ ■ ■ <g> Ut ® U ® E/i ® U ® • • • , 17 = -p • 1 > = ~P ■. ■ • ( L4 ) 





In this way, it is possible to use the VOA results in the principal regime in order to find form-factors in 
the disordered region. It remains only to take the limit J x — > J y limit, in order to obtain form- factors 
for the massless XXZ model. 

However, things are not quite so simple: while it is true that the VOA to the XYZ model par- 
allels that of the XXZ model, it does differ in one important respect. The explicit multiple- integral 
expressions for form-factors in the XXZ case are obtained by using a Bosonization technique - more 
precisely a free-field-representation of the quantum affine algebra U q (s{%). The problem is that such a 
free-field representation has not yet been found for the quantum elliptic algebra Aq lP (sl2) relevant to 
the XYZ model. The reason for this technical problem is ultimately linked to the absence of charge 
conservation around a vertex in the 8-vertex model associated with the XYZ chain. 

This problem has been considered before in the literature, and there are two ways to get around 
it. The first approach involves mapping the XYZ model to a solid-on-solid (SOS) model, for which a 
free-field realisation does exist (and for which the relevant algebraic structure is the elliptic quantum 
group of face type ^^(sb) [Ml ESSE])- This was the method developed and used by Lashkevich and 
Pugai to obtain expressions for both form-factors and correlation functions in the principal regime 
XYZ model |37l [38] . The method was extended to higher spin analogues of the XYZ model in [39 j . 
The second approach is applicable specifically to the massless XXZ model; the idea here is to derive 
and solve a difference equation (a deformed KZ-equation) for correlation functions from the analogous 
equation for the XYZ model [40J, or to construct a field realisation only after having already mapped 
to the disordered regime and taken the massless limit [27] . 

In this paper, we take the first approach, the main reason being that simplified expressions for 
the resulting XYZ two-particle form-factors mapped to the disordered regime are already present in 
the literature |41j . Our contribution is to take the appropriate massless XXZ scaling limit of these 
existing results, and to use them to compute the exact two-particle contribution to the longitudinal 
structure factor. 

In Section 2 of this paper, we describe the key components of the VOA to the XYZ model in 
the principal regime, the map to the disordered regime, and the limit from the existing disordered- 
regime XYZ form- factor results to our new massless XXZ form- factor expressions. In Section 3, we 



give the derivation of expression (3.4) for the two-spinon contribution to the longitudinal structure 
factor S zz (k,w). This was the key result quoted in the earlier paper |15| . In Section 4, we present a 
detailed quantitative analysis of the structure factor, an analytic derivation of the asymptotic threshold 
behaviour close to the upper and lower w limits of the two-particle (w, k) continuum, and an analysis 
of the region of the (w, k) band over which this threshold behaviour is a good practical approximation 
for different A values. We present some concluding remarks in Section 5. Finally, in Appendices A 
and B, we give the definitions and required properties of elliptic functions, and present an alternative 
derivation of the mapping of principal form factors to disordered ones. 
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2 From XYZ to Massless XXZ 

A general multiple-integral expression for principal regime XYZ form-factors can be constructed by 
following the approach of [371 EE] • The case of the form-factor of the operator a z is considered in detail 
in the paper [38j . and the author demonstrates a technique that enables him to obtain an expression 
for this form- factor which involves no integrals. This approach is extended to a x and a y form- factors 
by Lukyanov and Terras in fH]. Using the mapping mentioned in the previous section these authors 
present results directly in the disordered region of the XYZ model. In this section, we review these 
results and take the appropriate scaling limit to the massless XXZ model. This limit is different to 
the sine-Gordon limit considered in |41j . 



2.1 The XYZ model in the principal regime 

The XYZ Hamiltonian is derived from the 8-vertex model elliptic R matrix [33] given by 

( a{u) d(u)^ 



R(u) = p(u) 



b(u) c(u) 
c(u) b(u) 



\d(u) 



(2.1) 



i{u)j 



with 



a(u) = snh(A(l — u)), b(u) = snh(Au), c(u) = snh(A), d(u) = /csnh(A(l — u))snh(Au)snh(A). 

Here snh(tt) = — isn(iu) is the Jacobi elliptic function with modulus k. The definitions, relations 
between and required properties of all elliptic functions used in this paper can be found in Appendix 
|A} Let K, K' be the corresponding complete elliptic integrals given in Appendix|A| We use the variables 



„2r 



_ ttK' 

e k 



e in , C = x , 



_4irJf ■> 

from which definition it follows that A = K jr. We also define p = e k' and 5 = and sometimes 
use £ = r — 1 when connecting to the results of [H]. The principal regime is given by < x r < x < 
C < 1. We choose the scalar function p(u) as follows. 



p(u) 
(z;qi, • 



!- r/ 2 (* 4r ;* 4r )oo QM^x^eM^C 1 ) g(C 2 ) 

g(t 2 ) 



®Mx 2 (- 2 ) 

„2r\ 



(ry> 4 * ry* 4 ryi 2 V \ ( ry*2if* y ■ ' r y>4 rf^/V *\ ' 

-6,0- ,0- ^COV ? ? /CO 

OO 

wi,— ,rt m =0 

= {q;q)oo(z;q) 00 (q/z;q) 00 . 



This form of the R-matrix coincides with that used in [38] except for a minus sign in d(u). However, 
our notation differs slightly: most importantly, our p is not equal to the p of [38] ■ A full dictionary 
between our notation and that of both [38] and |41j is given in Table [!} 
For V = Cv + © Cv_ , we regard R(u) as a linear map on V <8> V by 



R(u)v ei <g> t> 



= 2 



£ , 1 ,£2=± 



R £ J £ J(u)v F 
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Table 1: A Dictionary of Notation 



This paper 


Reference [38] 


Reference [H] 


u 


u 




A 






r 


r 
2e 

7T 


l 


£ = r- 1 


r — 1 


e = _a_ 

1— ?7 


•n-A 

x = e 2A- 


x = e _e 


e -7r5/2 


2r 


2r 

p = X 


e -T<5(e+i) 


p = e k> 


— Z7T 

e «■ 


47T 



Let 14 (z = 0, 1 , . . , iV) denote A + 1 copies of V and regard Rij (u) as a linear map on Vn ® • • • <8> (8) Vq 
acting on the z-th and j-th tensor components as R{u) and on the other components trivially. We 
define the finite transfer matrix T(u) by 

T(u) = ti v Rno(u) R N - 10 (u) ■■■ R w (u). 



Then one can verify that in the infinite N limit the XYZ Hamiltonian (1.3) is obtained as 

JsnfA, k')cn(X, k') d , m . . , , 

H XY z = ' ' — hiT(u) +constant, 2.2 

zA ait m=o 

with k' = \/l — k 2 and 

J x = j(cn 2 (A,/fc') + /fcsn 2 (A,A; , )) , 

J y = J (cn 2 (A, fe') - fcsn 2 (A, A/)) , (2.3) 
J z = — Jdn(A, k). 



2.2 The vertex operator approach to the XYZ model in the principal regime 

The vertex operator approach to the infinite-lattice massive antiferromagnetic XXZ model was de- 
veloped in [22] and is described in detail in the book [20]. This approach was then extended to the 
principal regime XYZ model in |28j-|32j. The essence of the approach is to identify the transfer matrix, 
the space on which it acts, and local operators in terms of the representation theory of the underlying 
symmetry algebra. In the XXZ case, this algebra is the quantum affine algebra Uqisi^)] in the XYZ 
case it is the elliptic algebra Aq^sh) (beware that the (q,p) indicated in the name Aj.pC 5 ^) do not 
correspond directly to the notation of the current paper - in fact we have (q,p) = (—x,x 2r )). In this 
section and Appendix |Bj we give a brief outline of the key features of the approach relevant to the 
present work. We refer the interested reader to the original articles cited above and [39J for further 
details. 

The lattice transfer matrix of the infinite-size principal-regime XYZ model acts on the infinite 
tensor product space •••(g>F<8>^<8>^<S>V r (X>-- - with antiferromagnetic boundary conditions at plus 
and minus infinity. Let us choose to label the position of our spin-chain sites as • • • ,2,1,0, —1, —2, 
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Then the two antiferromagnetic boundary conditions we consider are labelled by j = or j = 1 and 
correspond to considering only those states which have the spin at site i, denote e(i), restricted to 
e^\i) = (— for \i\ » 0. More precisely, we introduce the space of states W^' (j = 0, 1) as a 
half of the infinite tensor space with the antiferromagnetic boundary condition j . Namely 

H ij) := Span c {• • • <g> v e{1) v e{0) \ e(i) = ±, e{i) = e {j) {i) for i > o} , (2.4) 

where v+ = , V- = • ^ ne starting point of the vertex operator approach is to identify the 

space of states with the level- 1 highest weight modules V(A?) U = 0> 1) of the algebra Aj^C 5 ^), 
where Aj denotes the fundamental weight of s^- Then, the total infinite tensor product space is 
identified with the tensor product 

jV) ■= n {j) n* ij) ~ End(^ (j) ). (2.5) 

The transfer matrix of the XYZ model and local operators are then identified in terms of certain 
vertex operators that act on the space (2.4). The relevant 'type I' vertex operators are maps involving 
both H^) anc l a finite-dimensional Aq tP {sl2) evaluation module V u = V®C[(, C -1 ]- They are -4.q,p(sb) 
homomorphisms of the form 

^-^(u) : n ij) — ► H {1 - j) ® V u , (2.6) 
whose components $>± (u) are defined by 

= ^~ j ' j) {u) ® v £ . 

e=± 

The transfer matrix of the XYZ model is identified with the map T{u) : -> defined by 

T(n) = ^"^(n) ® ^i'^iuf, (2.7) 

s=± 

where i denotes transpose. The XYZ Hamiltonian is then identified as 

Jsn(A, fc')cn(A, k') d 



2A du 



lnT(u) . (2.8) 



u=0 



Now we consider the 2x2 matrix E ee i that acts as E £)£ iv a = b a ^v e at site of the lattice. This 
local operator is realised as an operator on in the vertex operator approach as 

0(E ££ ,)^ = ^- j \u-l)^,-"\u) . (2.9) 



14=0 



Then the spin operator a x for example is realised as 0{a x )^' = 0(E^ + 0(E ^)^. 

The vacuum eigenvector of the transfer matrix T(w) is denoted by \v&c;pr)^ 6 J-^ and defined 

by 

T(«)|vac; pr) w) = |vac; pr)* 1- ^. (2.10) 
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This eigenvector^ has a very simple form in the vertex operator picture, and is constructed in terms 
of a certain grading operator H^> that acts on H^K More precisely, we define = — \p-\- 1, where 
p = Ao + Ai , and identify 

\va,c;pr)^ = r x 2H<,) , 

(Zti))3 

where we are regarding — End('Hw). Namely, T(«) acts on / £ End(7^ J ')) as 

T(u)/ = ^ $(W-i)( u ) o/o c&^'V), (2-11) 

see |2U| . The normalisation is defined by the 5I2 principal character 

zW = Tr w , ) (x 4 ^ ) ) = ? -^ r , (2.12) 

and is chosen such that w)(vac;pr|vac;pr)^^ = 1. Here the inner product of two elements f,g G 
End(7£Ci')) i s defined by (f,g) = Tr nU) (f o g). We denote a vector in = ftCi) $ by a ket 

vector in this section, but it should be understood that it is identified with an operator in End^^)) 
whenever one considers an action of the vertex operators on it. In what follows, we refer to arguments 
based on the identification T^' = End(H^) as the 'vertex operator picture'. 

The reader may at this point be thinking that the vertex operator approach is wholly algebraic and 
formal, but in fact <& e (u) and Z^ have a direct lattice interpretation in terms of the 8- vertex 

model: the operator is identified with Baxter's corner-transfer-matrix Hamiltonian, 5> e (u) with 
the half-transfer matrix, and Z^' with the partition function. In fact, it was Baxter's observation [33J 
that it it possible to express the partition function in terms of the corner-transfer-matrix Hamiltonian 
as in Equation (2.12), and the subsequent observation that this partition function was related to 
the sl2 principal character, that were the starting points for the development of the vertex operator 
approach. 

In order to construct other eigenstates of the operator (2.7) it is necessary to introduce a new 'type 
II' vertex operator ^f*^ 1 "^'^ (u), defined as the map 

^*d-i,i)( u ) . Vu ® n {j) — > n {1 - j \ 

with components XE^ 1 (it) specified by 

The full list of properties of both type I and type II vertex operators can be found in [38]. One 
property that we require in the current paper is the commutation relation of type I and type II vertex 
operators: 

$tt,l-j)( Ul ^*(l-jj)( U2 ) = t{ui _ U2 )**pM-i)( U2 )$(l-i,i)( Ul ) ) (2.13) 
where the function t(u) is given by 



01 (\-^pA 

r{u) = i^- ± (2.14) 

0i(| + f,^ 



1 We are using the term eigenvector loosely: |vac;pr)^ is a true eigenvector only of T 2 (u). 
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Our convention for theta functions is given in Appendix [XJ 
Let us consider a state denned by 

1 0! , 2 ; pr)g ea = m*Jj< l -ti (W 2 M Vti 1 '^ (t0i /tt) |vac; pr) ^ . 

It then follows immediately from (2.13), ( |2.7[ ) and ( |2.11 ) that we have 



r(u)|^,0 2 ;pr)g £2 = t(u - iOtMriu - i6 2 /ir)\9 1 , d 2 \pr)^g 



(i-i) 



(2.15) 



Hence, we may create a new eigenstate of T(u), i.e., an excited state with an eigenvalue t(u — 
i9i/7r)T(u — 162/7?), by acting on the vacuum with the type II vertex operators. More precisely, the 
type II vertex operators 'I'*^ 1 J (u) are identified with the creation operators of quasi-particle (spinon) 
excitations with rapidity 9 = —iixu and spin e. The 2n-spinon state (spinons are always excited in 
pairs) is given by 



n, 



~>2n 



U) 



(i8 2n /rr) ■ ■ ■ (i6i/ir) \vac; pr) ij) . 



(2.16) 



From (2.15), one finds that the eigenvalue of T(u) per spinon is r(u — iO/tt). Then from (2.8), one 



can deduce that the pseudomomentum k(9) and energy oj(9) of a spinon state in the principal XYZ 
model are 



Hence we obtain 



e ik{9) 

u>{9) 

k{9) 
u(9) 



Jsn(A, k')cn(X, k') d 
2A !hi 

21' e \ tt 



lnr(u — i9/ix) 



u=0 



am 



JI'snfA, k')cn(X, k') . (21 '9 , 
-dn I , kj 



A 



7T 



(2.17) 



(2.18) 



A 

Here we have introduced new complete elliptic integrals /, /' by 

x = e 1 

and denote by ki, k'j the corresponding moduli (that is, we now consider elliptic function with nome 
x as opposed to the original functions involved in the R-matrix which had nome x 2r ). The symbols 



sn(u, ki), dn(ti, ki), am(u, kj) denote Jacobi's elliptic function with modulus kj. In deriving (2.17) and 



(2.18), we have used the identity between elliptic functions of different nomes given in Equation (A. 7) 



Expressions (2.17) and (2.18) for spinon pseudomomentum and energy are consistent with the results 
of 



The 2n-spinon form factor of the local operator E ee i can now be expressed in the vertex operator 
picture as the following trace: 



00 



(vac;pr\E £)£l \9 

±tr H(J) ^O^tfl^H^fN^^/j) . . . yf-^^/n)) . (2.19) 
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The massive, antiferromagnetic XXZ model corresponds to the r — > oo (k — > 0) limit of the above 
picture. In this limit the above trace can be computed directly in terms of the free- field realisation 
of U q (sl2) [20] • However, there is no known free-field realisation for the general elliptic case. This 
problem was overcome in [37, 38J by mapping the 8-vertex model to the SOS model using Baxter's 
intertwiners. A free field realisation does exist for the SOS model [HI HSJ [Ml G5] and this was used to 



produced an integral expression for (2.19) which may be found in [38 



2.3 The map to the disordered regime 

Any regime of the XYZ model can be obtained from the principal regime by a suitable gauge trans- 
formation |33j . In this section, we construct such a transformation in terms of the following matrices 

^(i-O'^M- 1 -)- 

The adjoint action of these matrices on Pauli matrices is given by 

U (a x ,ay,a z )U^ = (aW*), U x {a x , a\ a z )U^ = {a\-cy\-a x ). (2.20) 
We consider the following gauge transformations. 

R(u) = (Ut ® U )R(u){Uq 1 ® u i l ) = {U Q ® U!)R(u)(U{ 1 ® U^ 1 ). (2.21) 



Note that a similar gauge transformation has been discussed in [27] . The difference is due to the shift 
A -> A - 2iK made in sec. 2.4 of [27] . 

Now define the infinite tensor product Gj (j = 0, 1) by 

Gj = ■••<& Ui-j (g> Uj (g> Ui-j (8) Uj (8) • • • . (2.22) 

where Uj acts at even sites of our infinite product space. Then it follows that with Hxyz given by 



(1.3), we have 

Tjdis /~i u s~t— 1 1 \ "V jdis x x , jdis^yv , jdis z z \ 

h xyz = GjHxYzGj =-jl^^ Jx (J i <T i+l + J V a i a i+l + J z ^i+l)' 



4 

where jf s = -J z , jf s = J x , jf s = -J y . 



With the parametrisation (2.3), we have \ J y \ < J x < —J z which corresponds to the principle regime. 
Hence, we have \J z ls \ < Jy ls < Jx ts which corresponds to the disordered regime [33]. 

In order to apply the gauge transformation to the space of states {I = 0, 1), let us devide Gj 
into two parts in the following way. 

Gj = Qj ® Qx-i 

with 

Qj = -- -®Uj% Ui-j (g> Uj, Gj = Uj (g> Ui-j g) Uj (8) ■ ■ ■ 
Here we assume the rightmost Uj of Gj acts on the 0-th site of our infinite product space. 



9 



By transforming by Qj, one finds the following two spaces. 

?4ia := S P an c{- ' • ® w ® W ® W <8> • • • <8> <8> Wj | io,ii, • • • e {0, 1} }, 
:= Span c {- • • <8 w\ <8 wi (8) ii>i <8 • • • <8 <8 u/j | jo, ji, • • • G {0, 1} }, 

where u>o,u>i denote the eigenvectors of cr 1 given by 
Namely we have 

^l = g j+ i-H^ (£ = 0,1). (2.23) 

Here and hereafter index j + £ should be understood in mod 2. We regard (j = 0, 1) as the 
spaces of states in the disordered regime. We also set H*JfJ = g 1+j+e n*^-^ and define the total 
spaced = GM^=H®e>H$. 

Accordingly, eigenstates of H^ Y z are obtained by acting with Gj on the eigenstates of Hxyz and 
have the same energy. Noting the duplication (2.23), we have a new vacuum vector |vac)^ in J 7 ^] 
expressed in two ways as 

|vac)^ = Gj + i\va,c;pr)( l ~^ 

with I = 0, 1. 



2.4 The massless XXZ limit 

The massless XXZ Hamiltonian is obtained by taking the limit x — >■ 1 (K — > +oo) of the disordered 
Hamiltonian Hxyz while keep r fixed. This corresponds to the following limits of the various elliptic 
parameters: 

2' ' ' + ' 2r 

Corresponding to this limit, we have 

Jx — > J, Jy —5" J COS ^ — ^ , J z — > — J. 

Defining A = cos (~), we then have 

which is the Hamiltonian of the massless XXZ model. 

In order to obtain the dispersion relation for the massless XXZ model, let us set 8+ = f3 and 

<5 



One should note pi = e~f . Then it follows from 12.171) and (12.181) that we have 



fc(0±) = ami— ,fc/)±^ 



7T 



JJ / sn(A 1 fcV(AjO /2J'/3 
w(0±) = r dn 

A V 7T 
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The above massless XXZ limit implies 

/^+oo, 

Defining k(/3) by 

k(/3) := 2arctan(e /3 ), 

we find the pseudomomentum k(8±) and energy w(0±) of a spinon of the massless XXZ model are 
given by 

limJfe(0 + ) = «(/3), limfc(0_) = k(/3) -tt 

limw(0±) = — r-r-r- = Ui?|sin/c(/3)|. 
cosh(p) 

Here vp denotes the Fermi velocity given by 

Jr /7r\ irJ a/1 — A 2 

WF = _ s in(- = — . 2.24 

2 Vr/ 2 arccos(A) v ' 

In deriving these limits we have made use of the conjugate modulus transformation for dn given by 



Equation (A. 5) and of the limits of elliptic functions given by Equations (A. 8) and (A. 9). 

For (5 real, the range of 2arctan(e' 3 ) is (0,7r), and so the + parametrisation gives us right-moving 
spinons occupying half the first Brillouin zone, and the - parametrisation gives us left-moving spinons 
occupying the other half. However, as we shall discuss in Section 3, spinons come only in pairs, and 
right-moving spinons are alone sufficient to span the complete Hilbert space of the quantum spin chain. 

2.4.1 The sine-Gordon Limit 

In [II], the sine-Gordon theory is discussed by taking a similar x — > 1 scaling limit of the XYZ 
model. It is interesting to compare our massless XXZ limit with this relativistic field theory limit. 
The approach of [H] involves shifting the rapidity in a different way: as 6 = $ — ? (we here use i? 
to indicate the parameter denoted by 9 in [U]). Then, with a slightly different normalisation of the 
Hamiltonian, the momentum p and the excitation energy e of the quantised soliton are given in [41j 
by 



# 4 (l _)- i^plj 



Then in the same scaling limit x — > 1 as above together with the limit that the lattice spacing e — > 0, 
the dispersion relation for massive relativistic particles is given in |41j as 

lim = M sinh <&, 

e 

£(■&) 

lim^^ = Mcoshtf, 
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where the mass M is given by 

7T 

M = lim . 

e 

This sine-Gordon limit is different to our massless XXZ limit and was used in [UJ in order to connect 
lattice and field theory operators. 



2.5 Form factors in the disordered regime 



In the same way that vacuum vectors in the disordered regime were constructed in Section 1.3, excited 
states of Hxyz are gi ven by 



|#i> • • • j #2 



in -Trfii with 



0,1. Hence a form factor of a local operator acting on the site for the disordered 



regime is given by 

^(v & c\E EjE ,\e 1: ■ ■ ■ ,e 2n ) { £,.., £2n = ^(vaciprlU^E^Uj+elOi,--- , 2n ;pr)C--W (2.25) 
The correspondingly gauge transformed Pauli operators acting at site of the lattice are given by 



A derivation of the formula (2.25) in the vertex operator picture is given in Appendix |B[ 

As we mentioned above, a general integral formula for the principal form factors appearing on 
the right-hand-side of Equation (2.25) can be found in |38j . Furthermore, this integral was there 

This same 



performed explicitly in the case when the local operator was a z - see (3.14)-(3.16) of 
method was then used to compute the integrals associated with the other Pauli operators in [41J and 
a summary of all cases can be found in Appendix A of [41J (where the results are presented directly 
in the disordered regime). Let us define the function 



f(zi,z 2 ) a A := 



1+1. 



2 



(2.26) 



where Fq and G(z,p) are defined in Appendix A of (UJ, and we now use the notation r = £ + 1. The 
results of Appendix A of |41j may be expressed succinctly as 



(?) 



(vac\<r c \0i,n2/± tT 
^(yac\ay\9x,e 2 )^ ± 
(^(vacl^l^,^)^ 



-m 



n T o 2 + ^±(-iy^f(e 1 + ^e 2 + n -)l% 



(2.27) 



*R)7(0i + ?, + j)% ± if (e, + > 2 + 



All other components, for example (vac\a z \6i, 6 2 )^P_, are zero. Note that e = ± labels on the 
form- factors (2.27) are inherited from spin labels in the principal regime but no longer have this 
interpretation in the disordered phase. Note also that the apparent ? shifts in the arguments of the 
functions / relative to |41| are again due to the fact that that our 6 and the corresponding symbol 
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in [41 j . which we here denote are related by 8 = $ — ?. We can remove these shifts by using the 
half-period property of theta functions given by Equation (A.6), from which it follows that 



7T 



Using this property leads to the following expressions for the two-spinon XYZ form-factors in the 
disordered regime: 



U) (v ac \a x \e 1 ,e 2 ) i i] T 



(-i)V(0i,02)H±/(0i,0 2 )K, 



(2.28) 



^Vck^A^l = (-i) j+1 /(^i^2)?;2T/(0i^ 2 )?;3- 

2.6 Massless XXZ form factors 

We now consider the disordered XYZ form factors in the x — > 1, r fixed, limit discussed in Section 



2.4 



For right-moving spinons, we identify 6 = (3 and take the p — > limit of the function /(/3i, /?2, c 
This limit is obtained from that of the theta functions and from the following limits: 

disinh((£ + l)t) sinh 2 ((l + ^)t) 



a,b 



KmG(p,p) = G(J3) := exp 



t sinh(^) sinh(2t) cosh(t) 



, limFo = 2(1 +r L )p s« . 



We find that only one of the six /(/3i,/3 2 )cd appearing in (2.28), namely /(/Si, ^2)12, is non-zero in 
this limit. As a result, the only non-zero, two-spinon massless XXZ form- factor is given by 



\2,3 



lim U) {yac\a z \l3 1 ,f3 2 } { J 1 ) 



-1V +1 £ 



ei,S2 



2 sin 



ft, _ ZL 
2i 4 



sm 



ft 
2; 



cos 



2i£ 



By an abuse of notation, we continue to use 
of this form factor. With /3i, 02 real, we have 



^ \vac I cr z I /3i, /S2)i^e 2 to refer to the massless XXZ limit 



2( l + r l)2 e -/ 4 (^) 



cosh(/3i) cosh(/3 2 )(cos(|) + cosh(ft^)) ' 



where the integral I^{z) is defined by 

dt sinh((£ + l)t (cosh(2i) cos(4tz) - 1) 



hi*) 



(2.29) 



(2.30) 



t sinh(£t) sinh(2t) cosh(t) 

Note that the expression (2.29) is symmetric with respect to exchange of (Si, ^2)) and independent of 
both e and j. 



3 The Longitudinal Structure Factor 



The longitudinal structure factor has already been defined by Equation (1.2). In order to be able to 
compute the form- factor expansion of this object we need to know the resolution of the identity in 
terms of a basis of states. We conjecture that 

I= ^Z] Yl f d ^"' I 'y^ 1 '"' >^ 2 ")K-,e2n ei.-iii^l'"- >#2»|- 

i=0,ln>0 {ZU) - J -°° Z7T J -°° Z7T 
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This conjecture is an analogue of the conjecture made in the antiferromagnetic regime of the XXZ 
model in [20]. Note however, that there is a slight but important difference to [20J. In considering XXZ 
single spinons in Section 2.2, we have characterised right-moving spinons as having momentum in the 



range (0, 7r) and left-moving spinons as having momentum in the range (— vr,0). In (3.1) however, we 



only include even spinon states consisting of right-moving spinons. The reason for this are two-fold: 
firstly, that right-left moving pairs are simply absent from the Bethe Ansatz states, and secondly that 
a right-right pair already spans the entire (0, 2tt) Brillouin zone - to include the left-left pair would 
be to double count. This point is discussed in some detail in |26j for the massive antiferromagnetic 
phase of the XXZ model. Independent numerical justification that we have made the correct choice 
of normalisation is given by the sum rule calculations in Section 4. 



In writing (1.2) we have not specified to which vacuum we are referring. However, the result is 
the same whether we choose to use |vac)(°), Ivac)^ 1 ) or the linear combinations |vac}± 
Ivac)^ 1 )) considered in |41j . For notational convenience, let us specify that |vac) = 



4(|vac>(°>± 



vac)(°). Then 



inserting the resolution (3.1) we have the 2-particle contribution 

/•oo 

■^rj, J OO 



s?{k,u) = -T / 

^ J — oo J — oo 



a pi ap 2 



where the spinon momentum and energy are defined by 

K(fi) = 2arctan(e /3 ), W{fi) 



cosh(/3) 

We then write both the j sum and w integral in terms of delta functions to give 

1 

4 



dfrdPi K k - K (A) - K (fc)) - W(pt) - W(f3 2 )) F(/3i,/3 2 ,£), (3.2) 



— oo J —oo 



where F(f3 1 ,(3 2 ,0 is defined by ^29 ). 

We now recall that if we have a suitably smooth function g : 
zeroes (x^ ' , x% P) at which detg'(x^ ,x\ r) ^ , and a function / : 



/oo roo 
/ dxidx 2 f(xi,x 2 )5(g(x 1 ,x 2 )) = ^ 
-oo J — oo 



/(4 



I 2 - 

l 2 - 
(j) 



I 2 with a finite number of 
I, then we have 



det g'(x\ 



To use this fact on (3.2), we make use of the determinant 



K'(lh) K'((3 2 ) 



VF(ta,nh((3 2 ) — tanh(/3i)) 
cosh(/3i) cosh(/?2) 



For each choice of k and to in the two-spinon band there is a unique (up to exchange of $\ and j3 2 ) 
pair (/3i,/?2) satisfying the combined conditions k = K(f3\) + K(f3 2 ) and oj = W{f3\) + W(f3 2 ). Thus, 



(3.2) becomes 



lcosh(/3i)cosh(/3 2 )F(/3i,/3 2 ,0 
2 



vf\ tanh(/?i) — tanh(/3 2 )| 

Let us denote the upper two-spinon energy threshold and and lower two-particle energy threshold 
by uJ2,u(k) and u 2t i(k). They are given by the following expressions: 

w 2 , n (fe) = 2vFsin(k/2), ui 2 j(k) = vf\ sin(fc)|, (3.3) 
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where the Fermi velocity vf is given by Equation (2.24). There is a useful identity 



w 2 u(^) — = vp\ tanh(/3i) — tanh(/?2)|. 



Using this identify, and substituting the expression (2.29) for the modulus squared of the form-factor, 
we arrive at the expression 



S z 2 z (k,u) 



@(uj 2 ,u{k) ~ w)0(w - W2,,(*i))(l + r 1 )^ 



yJul u {k)-U>* (cOs(|) + COsh(f) 

where /9 := ft - ft and fc = tf(ft) + A(ft), w = W(ft) + W(ft). 



(3.4) 



4 Results 

In this section, we offer quantitative results and plots for the longitudinal structure factor (LSF). 
Here, for convenience of comparison with previous results in the literature, we use a slightly different 
convention for the Hamiltonian (corresponding to the ones we used in [E]), writing it as 



(4.1) 



J > 0. We have a gauge equivalence H = OHxxzO^ , where O 



1 



<7 



where Sf = 2 

1 a 2 (g) • • • • The two-spinon part of the LSF is given by Equation (3.4) in which the parameter 
f3(k,u) =: 2irp(k,uj) is obtained from the constraint 



cosh(-7r p(k, w)) 



(4.2) 



4.1 Full Brillouin zone 



In Figures [T] and [2j we present plots of the two-spinon longitudinal structure factor for different 
values of anisotropy, starting at the XX limit and going up to the isotropic antiferromagnetic point 
A = 1 . All plots cover a full Brillouin zone and clearly show the continuum over which the two-spinon 
correlation is non- vanishing. This continuum in the k-u plane is located between the lower and upper 



boundaries (3.3). In fact, for zero magnetic field, the full LSF vanishes beneath the lower boundary, 
i.e., for uj < u)2,i(k). All 2, 4, 6, ... spinon states share the line to = u)2,i(k) as their lower boundary, and 
the LSF is strictly positive for all to > W2,l(k). Above the upper two-spinon boundary uj > ui2,u(k), 
the two-spinon contribution of course vanishes, but higher-spinon states can contribute. 

Starting at small A, the top left panel of Figure [T] clearly illustrates the fact that the LSF diverges 



at the upper boundary (we will quantify all threshold behaviour in Section 4.5), and tends to a 



constant at the lower one. This is easily understood [46| U7] by mapping to a system of free fermions 
using the Jordan- Wigner transformation, under which the S z operator becomes a fermionic density 
operator. Since the fermionic density operator in Fourier space takes the form of a convolution product 



of creation-annihilation operators, pk = V'fc+gV'fc; and since the ground state is a simple Fermi sea 
of the Jordan- Wigner fermions, all form factors of the S z operator are energy independent, and vanish 
for all but the two-spinon states. The LSF is thus simply a representation of the two-spinon density 
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Figure 1: Two-spinon part of the longitudinal structure factor of the infinite Heisenberg chain, for 
different values of the anisotropy parameter A (see also Figure [2]). For A — > 0, the correlation follows 
the density of states, and has a square root singularity at the upper threshold for all values of momenta. 

of states, which has a square-root divergence at the upper threshold (see, e.g., [H] for details) and is 
constant at the lower one. 

Turning the anisotropy up leads rapidly to the loss of the divergence at the upper threshold, 
starting from the region k ~ it, as can be seen in the top right panel of Figure [TJ The presence of 
a finite A is most importantly felt in the form factors: while the density of states still diverges, the 
form factors now decrease sufficiently rapidly to kill off this divergence. Once A has attained values of 
around 0.4, the remains of the divergence at the upper threshold have been completely erased for all 
values of momentum away from the zone boundary, and the lower threshold starts to feel the effects 
of the antiferromagnetic correlations, the peak at k = ir, uj = starting to develop. 

As A is increased further (see Figure [2j the cascade of correlation weight towards the low-energy 
sector continues until the isotropic limit is attained, at which point which most of the signal is con- 
centrated in the immediate vicinity of the lower threshold. 

Throughout this series of plots (which are presented in a uniform 00 and intensity scale for con- 
venience of comparison), the slow increase of the Fermi velocity and of the reach of the two-spinon 
continuum can be seen. 
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Figure 2: Two-spinon part of the longitudinal structure factor of the infinite Heisenberg chain, for 
different values of the anisotropy parameter A (see also Figure [2]) . Increasing the anisotropy shifts the 
weight progressively towards the lower boundary. The lower boundary becomes increasingly sharp as 
the A — > 1 limit is approached. 
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4.2 Fixed momentum 



In Figures [3]- [6j we provide a number of fixed-momentum cuts of the two-spinon part of the longitudinal 
structure factor. We organise the plots in each figure by increasing anisotropy at fixed momentum. 
Each individual plot also gives the A = and A = 1 curves as a reference; the effects of tuning the 
anisotropy continuously between these two values can thus be easily visualised. 

At k = 7r/4, as shown in Figure [3j the two-spinon continuum is quite narrow but its position 



is steadily increasing in energy as the anisotropy is turned up, since the Fermi velocity vf (2.24) is 



steadily increasing with increasing anisotropy. The upper-threshold singularity at A = disappears 



quickly as a function of A, being replaced by a square-root cusp (see the discussion in Section 4.5). 
The singularity at the lower threshold appears immediately but takes on a significant weight only for 
around A ~ 0.3 and above. The picture is very similar for the other momenta presented, namely 
k = tt/2 (Figure Q, k = 3vr/4 (Figure^ and vr (Figure [6]). 

4.3 Sum rules 

The full Hilbert space of the model in the zero magnetisation sector contains many more states than 
the simple two-spinon states we have considered. Two-spinon states in fact represent only a vanishingly 
small fraction of the total number of states when the system size goes to infinity. It is thus a remarkable 
fact that these simple states can carry a non- vanishing fraction of any correlation function. 

To quantify the importance of the two-spinon contribution to the longitudinal structure factor, we 
consider two sum rules. First of all, we use the integrated intensity 

obtained from the simple real space-time correlator (Sj, = j(t = 0)S*|(0)) = |. A less trivial sum 
rule comes from considering the integrated first frequency moment of the structure factor, giving the 
so-called f-sumrule (at fixed momentum) [49J, 

IfZ{k) = i ^(k,u } ) = --([[H,S z q ],St q }) = -2X x (l-co S k), (4.4) 

where 5| := ^= Y2j e%q ^ Sj 1S the Fourier transform of the local magnetisation operator and X x := 
{SjSj +1 ) is the ground state expectation value of the in-plane exchange term. The explicit value of 
the right-hand side of this identity can be obtained from the ground-state energy density eo |50j and 
its derivative using the Feynman-Hellman theorem, namely X x = 27(1 — A^)eo, with 



-•/(£ + !) 



eo = o sin 



7T 



00 , / tanht . 

dt [ 1 - z^urrr-r^ ) • ( 4 - 5 ) 

—00 \ 



tanh[(£ + l)t] 



The level of saturation of the two sum rules coming from two-spinon intermediate states is in- 
evitably anisotropy dependent. We provide the explicit values of the sum rule saturations coming 
from two-spinon contributions in Table [2] (for the f-sumrule, the saturation turns out to be exactly 
the same at all momenta). Two-spinon states carry the totality of the correlation at A = 0, a result 
which can be understood by considering the mapping to free fermions using the Jordan- Wigner trans- 
formation mentioned previously. Two-spinon states are the only intermediate states contributing to 
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Figure 3: Fixed momentum cuts at k = 7r/4 of the two-spinon part of the longitudinal structure factor 
of the infinite Heisenberg chain, for different values of the anisotropy parameter A. The A = and 
A = 1 limits are displayed in all plots for comparison. 
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Figure 4: Fixed momentum cuts at k = ir/2 of the two-spinon part of the longitudinal structure factor 
of the infinite Heisenberg chain, for different values of the anisotropy parameter A. The A = and 
A = 1 limits are displayed in all plots for comparison. 
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Figure 5: Fixed momentum cuts at k = 37r/4 of the two-spinon part of the longitudinal structure 
factor of the infinite Heisenberg chain, for different values of the anisotropy parameter A. The A = 
and A = 1 limits are displayed in all plots for comparison. 
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Figure 6: Fixed momentum cuts at k = tt of the two-spinon part of the longitudinal structure factor 
of the infinite Heisenberg chain, for different values of the anisotropy parameter A. The A = and 
A = 1 limits are displayed in all plots for comparison. 
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the longitudinal structure factor, and the sum rules are saturated to 100%. Our results are of course 
consistent with this fact. 

A more remarkable fact is that the subset of two-spinon states continues to play such a determi- 
nantal role in carrying the longitudinal structure factor even when the anisotropy has been turned on 
to significant values. As our results show, two-spinon states carry essentially all the correlation weight 
up to surprisingly large values of interactions A ~ 0.8, above which four, six, ... spinon states become 
harder to neglect. This is quite surprising since, reasoning again in the fermionic language obtained 
from the Jordan- Wigner transformation, the interaction should be able to create multiple particle-hole 
states quite easily, and arbitrarily complicated higher-spinon states should therefore participate in the 
correlation, leaving two-spinon states only negligible contribution. While true in the generic finite 
magnetic field case, it turns out that for zero magnetic field, the available particle phase-space shrinks 
to zero, and only the hole (spinon) part can disperse. The longitudinal structure factor thus possesses 
a finite two-spinon contribution, which is not true for example of the transverse structure factor. 
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0.7 


0.9637 


0.9578 


0.2 


0.9986 


0.9984 




0.8 


0.9406 


0.9314 


0.3 


0.9964 


0.9959 




0.9 


0.8980 


0.8844 


0.4 


0.9927 


0.9917 




0.99 


0.7918 


0.7748 


0.5 


0.9869 


0.9849 




0.999 


0.7494 


0.7331 



Table 2: Sum rule saturations as a function of anisotropy: two-spinon contribution to the integrated 
intensity l zz (4.3) and first frequency moment I[ z (4.4). 



4.4 Correspondence with finite size results 

The longitudinal structure factor can also be obtained at finite size using explicit summations over 
intermediate states, as performed in [51|. I52j . At finite size (i.e., when the chain is made of N sites; we 
take N even to ensure a non-degenerate ground state at zero field), two-spinon states in zero field are 
those eigenstates obtained using N/2 — 1 real rapidities and a single negative-parity one-string (those 
readers unfamiliar with this nomenclature are referred to the explanations in |52j). The total number 
of such states then corresponds to the number of ways of choosing two hole quantum numbers from 
N/2 + 1 available ones, and thus equals N(N + 2)/8. We can thus, at a fixed size N, sum over the 
contributions to (1.2) coming from these two-spinon states, and compare with our infinite-size result. 

Figure [7] presents such a comparison, done at a representative value of anisotropy A = 0.7 and 
two values of momentum, k = tt/2 and k = ir. Other values of anisotropy and momenta give quali- 
tatively similar plots. Besides the thermodynamic limit curve obtained from plotting the two-spinon 
contribution we have obtained, we also present the equivalent curves for three different system sizes, 
N = 256, 512 and 1024. The finite-size results must be smoothed with a gaussian, since the correlation 
function is then a sum of delta-function peaks split in energy by the mean energy level spacing. This 
smoothing can be sharpened at increasing system size, and this increasing sharpness can be clearly 
seen, e.g., at the lower boundary of the k = ir/2 plot. It is clear that the finite-size curves tend to the 
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Figure 7: Comparison between the two-spinon longitudinal structure factor for A = 0.7 and momen- 
tum k = 7r/2 (left) and k = tt (right), in the thermodynamic limit (solid curve) and at finite size. 
The finite size results are obtained by summing exactly over the set of N(N + 2)/8 two-spinon states 
at chain length N, and are plotted after gaussian smoothing of the delta function contributions. See 
main text for more details. 



thermodynamic one as system size increases; the inset of each plot offers a magnified view of a selected 
region away from the singular thresholds, the quantitative differences between the highest system size 
considered and the infinite size limit being of the order of a percent. 
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0.9664 


0.9636 


extrap 
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0.957(4) 


0.957(4) 


oo 


0.9637 


0.9578 


0.9578 



Table 3: Sum rule contributions obtained from the finite-size calculation of the two-spinon contribu- 
tions to the longitudinal structure factor at the representative value of anisotropy A = 0.7, for system 
sizes up to N = 1024 and extrapolated (see text), as compared to the analytical result at infinite 
system size. See main text for more details. 

In Table [3j we provide quantitative results for the sum rule contributions for the representative 
value of anisotropy A = 0.7. Similarly to Table [2j we provide both the integrated intensity sum rule 



contribution (4.3) as well as the f-sumrule one ( |4.4[ ), coming from two-spinon states. The results at 
finite size are clearly seen to tend to their infinite-size limit. For completeness, we have extrapolated 
the finite-size results using data at N = 256,512,768 and 1024, fitting with a polynomial in 1/^/N. 
Within the accuracy of the extrapolation, these results coincide with the ones obtained from the 
analytical form. 
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4.5 Threshold behaviour 



The investigation of the precise form of the longitudinal structure factor near the upper and lower 
boundaries of the two-spinon continuum, which our results render possible, is of considerable theoreti- 
cal interest in view of recent developments in the general phenomenology of one-dimensional quantum 
liquids |53[ 1541 [55] coming from efforts to calculate dynamical correlations away from the low-energy 
limit [561 E3 EHl EHl [6Ql EU |62] . For the special case of the longitudinal structure factor of the zero-field 
gapless XXZ antiferromagnet we are considering, which is equivalent to the density-density correlator 
of spinless interacting fermions obtained via a Jordan- Wigner transformation, the field theory predic- 
tions of |61[ [62] yield a singular LSF in the vicinity of the lower /upper thresholds of the two-spinon 
continuum. In the vicinity of the upper threshold, for < A < 1, this singularity is shown to be of the 
form S zz (k,to) rrt Wto~2 u(k) — to. For the lower threshold, the power- law becomes anisotropy- 

dependent, S zz (k,uj) 7^ (to — t02 t i(k)) ~( 1 in which K is the Luttinger parameter, which for 

the zero field XXZ model takes the value K = |(1 — acos In terms of the parameter £, the 

Luttinger parameter becomes K = ^(1 + |). Our discussion here has two aims: firstly, to reproduce 
and possibly refine the determination of this threshold behaviour; secondly, to quantify its region of 
validity, which is very difficult to achieve within nonlinear Luttinger liquid theory. 

Let us thus consider evaluating the two-spinon part of the longitudinal structure factor in the 



vicinity of the excitation thresholds, starting from our exact representation (3.4). We consider the 



explicit evaluation of the fundamental integral (2.30) in various limits. For convenience, we rewrite it 

as 



h(p)=lf ) (p)-lf\ P ) 1 



where 



r m. . f°° dt sinh(£ + l)t sinhi A ,„ „. 

€ y ' Jo t sinh^t sinhtcosh 2 i v ' 

Throughout the discussion below, unless we specifically mention otherwise, we consider the generic 
case £ = 0(1). 

4.5.1 The structure factor near the upper threshold. 

In order to obtain the structure factor near the upper boundary of the two-spinon continuum, i.e., 
to — > u)2,u(k), we must consider the limit p — > of the fundamental integral. Let us thus look at the 



two integrals (4.6) and (4.7) in turn 



a) I^(p): For p = 0, the integrand of 1^ is regular at t — > 0, but the integral diverges logarithmically 

as t — > 00. We thus expect a log divergence as a function of p, ii -(p) — > — chip + d where c > and 
c,d = 0(1). We can in fact immediately predict the value of the coefficient c by looking at the ratio 
of hyperbolic functions, which tends to 2 at t — > 00, so c = 2. This is easily proved by rewriting 

(defining f^\t) := ^^7^ cosh 2 ^ usm S t ne cosine integral Ci as 

4 1] ( P ) = 4 n) o#i) + 4 12 \p\h) - 2Ci(4 P t 1 ), 
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in which t\ is an arbitrary real number and 

4 n \p\h) •= r j fV(t) cos Apt, 4 12 Hp\h) ■= r - 2)cos4pt, 

and we have used the identity 




— cos Apt = -C\{Apt\). 



Let us choose t\ = 0(1). We then have 

l4 11) (p)l</ 1 f/f ) (*) = o(i). 

We thus have explicit convergence of the first integral, lim^o 1^ (p\ii) = 0(1). 
i~( 12 ) also converges explicitly for generic ^ = 0(1), 

l^ 11 ^/ )! ^ /°° - 2 )l = o(i). 

The only problematic terms as p — > are therefore relegated to the cosine integral, which can be 
rewritten to separate out the singular p dependence, 

ft A+ 

Ci(4pt) = C + In p + ln(4t) - 2 / — sin 2 2/ji. 
We thus obtain the partial result 

-21np + 0(l). 

6) lj?\p): The integrand of I 1 - 2 ) vanishes sufficiently rapidly at t — > and i — > oo, so this integral 
yields a contribution of order p 2 — > 0. Thus no nontrivial contribution to the structure factor comes 
from this integral in the limit considered. 



This means that overall, we have 

l c (p) -21np + 0(l). (4.8) 



We thus have (using p ~ y/w2,u{k) — ui from (4.2)) 

S ?^ u ) m){^ k 2 r 1/2 ^2Ak)-u (4.9) 

in which f u {C) is a momentum-independent function of anisotropy. The exponent we obtain confirms 
the field theory predictions |61j for the anisotropy-independent square-root cusp at the threshold (for 
< A < 1). Our results allow us additionally to extract a strongly momentum-dependent prefactor, 
which greatly enhances the spectral weight around the zone boundaries at k = 0,2n, as is also 
noticeable in the figures. 
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For the A — > limit (so £ — > 1), we have to take the limit more carefully, since the cosh +cos j 
in the denominator of the structure factor now vanishes when p — > 0. Overall, in this case one rather 
obtains a square-root divergence, 

(sin h\-l/2 

Sl z (k,co) jrt /u(l) ; 2 = (4.10) 

as expected, since in this non- inter acting case the structure factor simply follows the density of two- 
spinon states, all two-spinon form factors being energy independent and equal to each other. 

4.5.2 The structure factor near the lower threshold 

To evaluate the structure factor near the lower threshold of the two-spinon continuum, i.e., for u> —> 
^2,i{k), we need to consider the limit p — > oo of the fundamental integral. Here, we again consider 



£ = 0(1), and split the integrals precisely as before, using (4.6), (|4.7[ ). Again choosing ti = 0(1), we 
can see that 1^ (p) is still bounded by a constant, and so is If%>). In fact, since the integrands 
oscillate rapidly, this constant is zero. Moreover, the cosine integral evaluated at infinity also vanishes, 
so we have ii 1 ^ (p) > 0. 

(2) 

For 1^ (p), we start by writing 



/f ) (p) = /f 1) (p) + /f 2) (p), (4.11) 

where we have defined 



if ) (P ) = r if \ P) = r d ± sm22pt 2 . (4.i2) 

« yH> J t sinh2i ' « y " J t tanh£tcosh 2 i 1 ; 



We have 



/f \ p ) =lncosh7rp~vrp + 0(l). (4.13) 
We can evaluate 1^ (p) for large p by splitting it up, 

if V) = r g sin2 ^ 2 + r * ™ 22p \ (4.14) 



/o i tanh£icosh t Jt 2 t tanh£fcosh t 
Let us choose t<i = 1/y/p (any power between and 1 would do). For the first integral, we can write 

Wp dt sm 2 2pt 2p„ n , 2 ,. /"V? , sin 2 t ,„ ir . 

1 <l + 0(p- 2 )) dt—^. (4.15) 



o 



t tanh£tcosh 2 i £ Jo * 2 



We also have 

r2 v^ , sin 2 t vr 1 1 



In the second integral (from £2 = 1/y^ to oo), sin 2 2pt rapidly oscillates and we can thus replace it by 
1/2 when taking the limit p — > oo. This yields 

dt sin 2 2pt . u .l f°° dt 
~-T x = (l + O^ 1 )); 1 



'i/Vp * tanh£icosh i ' 2 Ji/^p t tanh£tcosh t 

(1 + 0(0) 



dt f 1 1 ^ Jf> 



2j X j^p t \tanh£tcosh t £t ) 2£ 



f + 0(1). 
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Adding up, we thus get 



r(2) > * / 1 



This means that the overall behaviour of the fundamental integral is 



We thus find (using p ~ ^ ln( ^_^ ^ ) from (4.2)) that the structure factor behaves as 



sin ifer^W) (sin |)~f 



SZ*(k,u) ^tfi(0 \, 2 A , (4-18) 

[U - LJ 2 ,l{k)\ 2K ^ 

where fi(£) is again a momentum-independent function of anisotropy. This result assumes that k ^ ir; 

i i 

at the antiferromagnetic point k = ir, the behaviour becomes 1/u) « due to the vanishing of U2 i- 
For the A — > limit (so £ — > 1), we thus get 

as expected, since the two-spinon density of states is simply a constant in this region of the continuum, 
and so are the form factors. 

Our analytical form for the two-spinon part of the structure factor has thus allowed us to reobtain 
the threshold exponents predicted from field theory, and to complement the threshold behaviour of 
the longitudinal structure factor with momentum-dependent prefactors hard to access within that 
method. 

One final comment here concerns the potential effect of higher-spinon states on the threshold 
behaviour. For the generic < A < 1 case, the obtained threshold exponents would remain unchanged 
upon the addition of these contributions, since the power series in the energy distance 5oj to the 
singularity cannot contain any logarithmic terms which could re-exponentiate into a different power 
law. The prefactor might however be corrected. 

4.5.3 Region of validity of threshold behaviour 

We now compare the behaviour at the lower and upper thresholds with the numerical evaluation of 



expression (3.4) in order to see over what range of frequencies these threshold formulae are valid. 
We do this in two different ways. The first way consists in evaluating the ratio S\ z I S?fi r at fixed 
momentum as a function of u, in which Sf? represents the relevant threshold behaviour in equations 



4.9|4.18 ) (the prefactors f u ,l{0 being obtained numerically directly from the exact representation 



3.4[)), and to find the region of u near the singularity for which this ratio remains one within the 



required accuracy. The second way consists in actually fitting a plot of the exact expression (3.4) with 
the expected threshold power law over a finite but small frequency region near the singularity, and to 
then check over which interval in frequency this fit remains consistent. The latter method emulates 
the kind of fitting one might do starting from approximate ab-initio numerical data for the structure 
factor, and gives an overestimate of the region of validity. The two methods give results consistent 
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with each other when the region of validity is at least of a few percent of the available continuum. 
Otherwise the stricter first way gives a much smaller region of validity. 

Tables [4] , [5] and [6] summarise where the difference between the threshold formula and the numerics 
becomes more than 1% and 10%, as a percentage of the numerical result, using the first criterion. 
The entries in the final four columns of these tables show the approximate percentage of the to range 
(u u — oji) that the threshold behaviour is valid for (within 1% and 10%). For example, for A = 0.5 
and k = 0.125 (the first row in Table [5]), the formula for the upper threshold is within 1% of the result 
for ~ 7% of the oj range, and within 10% of the result for ~ 47% of the oj range. 

It is immediately obvious that the region of validity of the lower and upper threshold behaviours 
depends strongly on the anisotropy. For low anisotropy, the upper threshold is not fitted well, the 
lower one being better described. For intermediate anisotropy the fitting is very reasonable, and covers 
a substantial range (over half) of the continuum. At high anisotropy, the lower threshold is not well 
fitted, whereas the upper one is rather well approximated. One point to notice is that the range of 
correspondence between the lower threshold behaviour and the exact structure factor for momenta 
at or near ir is very narrow: in this case, subleading terms correcting the theshold behaviour should 
not be neglected (the very low percentages presented in some entries in the tables should however 
be considered as indicative only, in view of numerical difficulties in evaluating the structure factor 
in the immediate vicinity of thresholds). On the other hand, the less sensitive second way of fitting 
gives acceptable fits over a wider range of frequencies. Plots for the structure factor accompanied by 
the threshold fits using this second way are given for these three values of anisotropy and momentum 
values in Figures [8 10 



fc/(27r) 


Wl 


Wu 


Lower 
< 1% 


Lower 
< 10% 


Upper 

< 1% 


Upper 
< 10% 


0.125 


0.751 


0.813 


~ 1% 


~ 12% 


< 0.1% 


~ 0.6% 


0.25 


1.06 


1.50 


~ 1.2% 


~ 14% 


< 0.1% 


~ 0.5% 


0.375 


0.751 


1.96 


~ 1.8% 


~ 20% 


< 0.1% 


~ 0.4% 


0.5 





2.13 


< 0.1% 


~ 0.9% 


< 0.1% 


~ 0.3% 



Table 4: Approximate validity of threshold results for anisotropy A = 0.1. These results are based on 
the first way of fitting discussed in the text. The low percentages are very approximate. 



jfe/(27r) 


Wl 


Wu 


Lower 
< 1% 


Lower 
< 10% 


Upper 
< 1% 


Upper 
< 10% 


0.125 


0.919 


0.994 


~ 16% 


~ 38% 


-7% 


~ 47% 


0.25 


1.30 


1.84 


~ 16% 


~ 38% 


~ 9% 


-52% 


0.375 


0.919 


2.40 


~ 11% 


~ 35% 


~ 13% 


-57% 


0.5 





2.60 


< 0.1% 


~ 0.1% 


~ 14% 


~ 50% 



Table 5: Approximate validity of threshold results for anisotropy A = 0.5. These results are based on 
the first way of fitting discussed in the text. The low percentages are very approximate. 
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fv I I Z7T J 


Wl 


w u 


Lower 
^ 1 % 


Lower 

<~. lU/0 


Upper 

<^ J. /o 


Upper 

<- -LU/O 


0.125 


1.07 


1.16 


< 0.1% 


~ 0.3% 


~ 1.8% 


~ 16% 


0.25 


1.52 


2.15 


< 0.1% 


~ 0.3% 


~ 1.8% 


~ 16% 


0.375 


1.07 


2.81 


< 0.1% 


~ 0.25% 


~ 1.7% 


~ 14% 


0.5 





3.04 


< 0.1% 


< 0.1% 


~ 1.2% 


~ 11% 



Table 6: Approximate validity of threshold results for anisotropy A = 0.9. These results are based on 
the first way of fitting discussed in the text. The low percentages are very approximate. 

5 Conclusion 

In summary, we have presented an analytical expression for the two-spinon contribution to the longi- 
tudinal structure factor of the XX Z chain in the gapless antiferromagnetic regime < A < 1 for zero 
field in the infinite size limit at zero temperature. Our results extend to this region previous results 
for the isotropic or gapped antiferromagnet. 

The question of the transverse structure factor remains difficult for the methods presented here. 
In the basis we are using, all transverse spin operator form factors vanish upon taking the gapless 
limit, and this points to the need for a resummation over states including macroscopic numbers of 
spinons, something which goes beyond current capabilities. Extending these results to the case of a 
finite magnetic field faces similar issues; the restriction to zero temperature is even more severe. We 
leave these questions open for the moment. 

Another important (but now feasible) extension to our work would be to consider higher-spinon 
contributions to the longitudinal structure factor. This was actually performed for the isotropic 
XXX antiferromagnet in the recent past |24j : one can expect that such a calculation would yield 
the longitudinal structure factor to around 1% accuracy for any value of anisotropy in the gapless 
antiferromagnetic regime in zero field. It would also allow us to further refine the determination of 
the threshold behaviour and of its limit of applicability. We will investigate these and other issues in 
the future. 
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(a) k/2-K = 0.125 



(b) k/2-K = 0.25 





(c) k/2n = 0.375 



(d) k/2n = 0.5 



Figure 8: Threshold behaviour for A = 0.1, the solid line is the numerical evaluation of S^, while the 
dashed and dotted lines indicate the upper and lower threshold behaviour respectively, fitted using 
the second method (see main text). 



A Elliptic Functions 

In this appendix, we collect together the notational conventions and identities for the various elliptic 
functions that we use in the paper. Many further properties of elliptic functions can be found for 
example in [63j . We firstly make use of the Jacobi elliptic functions defined by 



sn(u) = sin(am(n)), cn(u) = cos(am(n)), dn(u) = \J\ — k 2 sin 2 (am(ii)), 



(A.l) 
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(a) k/2n = 0.125 



(b) k/2n = 0.25 





(c) k/2n = 0.375 



(d) k/2n = 0.5 



Figure 9: Threshold behaviour for A = 0.5, the solid line is the numerical evaluation of S^, while the 
dashed and dotted lines indicate the upper and lower threshold behaviour respectively, fitted using 
the second method (see main text). 



where am(u) is the amplitude function defined in terms of the elliptic modulus k by 

»am(u) ^ 



yl — k 2 sin 2 (x) 



(A.2) 



Defining the conjugate modulus by kl = \f\ — k 2 , the complete elliptic integrals K, and K' are defined 
in terms of the same modulus k by the integrals: 



K 



dx 



a/1 — k 2 sin 2 (x) 



K' 



dx 



o \J\ - k' 2 sin 2 (x) 



(A.3) 
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1 .5 2 .0 2 .5 ^ 0.5 1.0 1.5 2.0 2.5 3fl 



(c) k/2-K = 0.375 (d) k/2ir = 0.5 

Figure 10: Threshold behaviour for A = 0.9, the solid line is the numerical evaluation of 5| z , while 
the dashed and dotted lines indicate the upper and lower threshold behaviour respectively, fitted using 
the second method (see main text). 

We sometimes make the k dependence of these various functions explicit by writing them as am(u, k), 
sn(u, k), cn(it, k) and dn(u, k). We also define the following functions 

snh(u) = — i sn(iii), cnh(u) = cn(iu), dnh(n) = dn(iu). (A. 4) 

The other type of elliptic function we use are theta functions. These are defined in terms of a 
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parameter q called the elliptic nome by 

oo 

$ x {u,q) = 2q 1 /\q 2 ;q 2 ) ooS m(Tru)Y[(l-2q 2n co S (2Tru) + q 4n ) 

n=l 

oo 

Mu,q) = 2q 1 / 4 (q 2 ;q 2 ) 00 cos(iru) Y[(l + 2q 2n cos(2iru) + q 4n ) 

n=l 

oo 

Mu,q) = (q 2 ;q 2 UH(l + 2q 2n - 1 C o S (2iru) + q 4n - 2 ) 

n=l 

oo 

Mu,q) = (q 2 ;q 2 )ooH(l-2q 2n - 1 co S (27ru) + q 4n - 2 ) 

n=l 
k' 

If we identify the elliptic nome as q = e n K , then the Jacobi elliptic functions and theta functions 
are related by 

sn(u) = — f^f — emu) = —j= — — r, dniu) = Vk' — — (-. 
In this paper, we use various identities that we now list. 

Identity 1 - the conjugate modulus transformation for Jacobi elliptic functions [63 : 

,. .sn(u, k') ,. . 1 .. .dn(«, k') rN 

sn(iu, k)=i — - — — , cn(tu,k)=i — - — — , dn(iu,k)=i — - — — . (A.5) 
cn(u,K) cn[u,K) cn(u,k) 

Identity 2 - the half-period property of theta functions |63j : 

M^± ~,q) = i+iy-e-^e^M^q), 

where q = e inT and 1 := 4, 2 := 3, 3 := 2, 4 := 1, 9l = g 4 = 1, g 2 = g 3 = 0. (A.6) 



* 7 ~ — TV =sn +icn ,ki). (A.7) 



Identity 3 - 



Identity 4 - the limiting behaviour of elliptic functions (which follows straight from the above defi- 
nitions): 

am(«, fc = 1) = 2 arctan(e n ) - — . (A. 8) 

am(u, k = 0) = u, sn(u, k = 0) = sin(it), cn(u, /c = 0) = cos(u), dn(u, A; = 0) = 1. (A. 9) 



B Derivation of (2.25) in the Vertex Operator Picture 



A key observation is an identification of the type I vertex operators with the half-transfer matrices on 



the lattice. Then applying the gauge transformations (2.21) to each R matrix constituting the half- 
transfer matrix, one can reach the following definition of the type I vertex operators in the disordered 
regime. 

: = E (Uj^Gw+izf^WTli ■■ Qi+fliw -+ g j+ t+iH®. (b.io) 

e'=± 
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Note that from (2.23) these are linear operators on Accordingly, one can realise a local operator 

as an operator on = Gj+tH^^ by 



0>{E Ee ,)f :=*^(«-l)*^(u 



0) 



u=0 



Then it follows from (2.9) and (B.10) that we have the gauge transformation of the spin operators as 

o'(a z p = (-y+% +e o(ay)^gj^. (b.h) 

This is consistent to the transformation of a a as a 2x2 matrix: AdC/" 1 : a x ,a y , a z \-¥ (— l) J cr z , a x , (— \y a y . 

In the vertex operator picture, the gauge transformation of the vacuum vectors |vac;pr)^' = 
x 2HW /(Zti)) 1 / 2 g S End(H^) follows from the fact that x 2hU) originates as the product of two 
corner transfer matrices. Hence the vacuum vectors as well as other physical 2n-spinon excited states 
should have the same transformation property as the corner transfer matrix. Therefore we identify 
new vacuum vectors as 



|vac)^ = Gj + £ , -.r 



•(^(l-/))l/2" 



2H( l - l ) r -l 



0,1) 



(B.12) 



in J 7 ^] = End(?{j^). Similarly to the type I case, the type II vertex operators Vl/ e {id/ir) in the 

(i) 

principal regime are mapped to linear operators on H^ is by 



(i-jj)/ 



One should note that the new type II vertex operators ^*^p(9) commute with the new type I vertex 
operators <&~1 (u) in pairs. Namely, 



Hence we obtain the following identification of the disordered 2n-spinon states 



v(j) 



7 l5 (7 2 • • ■t72n/ei,ea,-.sa 



»(3) 



7T 



rX 2ffC1 -^7 + V (B.13) 



Combining (B.ll), (B.12) and (B.13), we then obtain the form factor of the spin operator a a in the 
disordered regime as 



(vac|<7 a |0i, ■ ■ ■ , 02n)i{,- ,e 2n 



= ^(vaqprlAdCZ-V^)^!,-- - , «2„;pr)g^ eaB 
with £ = 0, 1. 



(B.14) 
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